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Quantum Search with Multiple Walk Steps per Oracle Query 


Thomas G. Won^ and Andris Ambaini^ 

Faculty of Computing, University of Latvia, Raipa bulv. 19, Riga, LV-1586, Latvia 

We identify a key difference between quantum search by discrete- and continuous-time quantum 
walks: a discrete-time walk typically performs one walk step per oracle query, whereas a continuous¬ 
time walk can effectively perform multiple walk steps per query while only counting query time. As 
a result, we show that continuous-time quantum walks can outperform their discrete-time counter¬ 
parts, even though both achieve quadratic speedups over their corresponding classical random walks. 
To provide greater equity, we allow the discrete-time quantum walk to also take multiple walk steps 
per oracle query while only counting queries. Then it matches the continuous-time algorithm’s run¬ 
time, but such that it is a cubic speedup over its corresponding classical random walk. This yields 
the first example of a greater-than-quadratic speedup for quantum search over its corresponding 
classical random walk. 


I. INTRODUCTION 

Quantum walks are the quantum analogues of classical 
random walks mm, and they have been the subject of 
much investigation for their algorithmic role in search [3] , 
element distinctness [3] , triangle finding [S] , and evaluat¬ 
ing NAND trees [6]. Formulated on a graph, a quantum 
particle walks locally from one vertex of the graph to 
adjacent vertices in superposition. 

As with classical Markov chains, quantum walks can 
evolve in discrete or continuous time. But these two 
formulations differ in the required number of degrees of 
freedom. With the N vertices of a graph labeling compu¬ 
tational basis states of an A^-dimensional Hilbert space, 
continuous-time quantum walks are well-defined on these 
vertices. Discrete-time quantum walks, however, require 
additional “coin” or spin degrees of freedom in order to 
evolve non-trivially laiH]. This leads to some algorith¬ 
mic differences between the two approaches 13 HQ], and 
much work has been done to show the relationship be¬ 
tween them 

In this paper, we identify a key difference between 
discrete- and continuous-time quantum walks as they are 
typically used in solving spatial search 113: where the 
goal is to find a “marked” vertex in a graph by querying 
an oracle. The oracle is considered to be an expensive 
“black box” that we want to use as little as possible, 
and other operations are “cheap.” Following this stan¬ 
dard for oracular problems, we compare algorithms by 
their oracle query complexity. As we will show in the 
following sections, the usual discrete-time quantum walk 
search algorithm takes one walk step per oracle query, 
while the standard dehnition of the continuous-time al¬ 
gorithm effectively allows multiple walk steps per oracle 
query. 

We show how this difference affects search on the “sim¬ 
plex of complete graphs,” which was first introduced in 
m, and an example of which is shown in Fig. In the 
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FIG. 1. A 5-simplex with each vertex replaced by the com¬ 
plete graph of 5 vertices. One complete graph is fully marked, 
indicated by double circles. Identically evolving vertices are 
identically colored and labeled. 


graph, we have M -|- 1 complete graphs of M vertices, 
arranged so that so that the vertices of a complete graph 
are each connected to different complete graphs. Here, we 
search for one fully marked complete graph, as indicated 
by the red a vertices with double circles in Fig. [I] We 
also have the blue b vertices that are one away from (ie., 
adjacent to) the a vertices, and the white c vertices that 
are two away from the a vertices. Classically, a random 
walker is expected to transition from one complete graph 
to another once every 0(M) steps, and it must make 
0(M) such transitions, on average, to hnd the marked 
complete graph, resulting in 0(M^) total steps. 

Next we analyze search on this by typical discrete-time 
and continuous-time quantum walks, highlighting that 
the number of walk steps per oracle query are different 
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and must be accounted for in showing quadratic speedups 
over classical. Then we adjust the discrete-time quantum 
walk so that it can also take multiple walk steps per oracle 
query. This yields a cubic speedup over its corresponding 
classical walk with the same number of walk steps per 
oracle query, and is the first example of a greater-than- 
quadratic speedup. 


II. DISCRETE-TIME QUANTUM WALKS 


We begin with the typical discrete-time quantum walk 
search algorithm The N = M{M + 1) vertices 
of the graph label computational basis states of an N- 
dimensional “vertex” Hilbert space, and the M directions 
that a particle can move from each vertex spans an ad¬ 
ditional M-dimensional “coin” Hilbert space. Together, 
C^(8)C^ is the Hilbert space of the system. Let |s„) and 
I Sc) be uniform superpositions over the vertex and coin 
spaces, respectively: 


N 


M 




y/M 


i=l 


Then the system \tfj) begins in 

IV'o) = |St>) ® |Sc), 


which is an equal superposition over both the vertex and 
coin spaces. The discrete-time quantum walk is obtained 
by repeated applications of 


Uo = S-{In<^ Co), 

where Co is the “Grover diffusion” coin [3] 
Co — 2|sc)(sc| Iai, 


and S is the “flip-flop” shift [10] that causes the particle 
to hop and then turn around, e.g., S'(|f) <8) \i ^ j)) = 
\j) (g) \j i). 

With this choice of initial state |^o) and evolution Uq, 
Fig.[T] shows that there are only three types of vertices, 
which we indicate by identical colors and labels, each 
with two types of directions. In particular, the a ver¬ 
tices can either point towards each other or towards the 
b vertices, the b vertices can either point towards the a 
vertices or towards the c vertices, and the c vertices can 
either point towards the b vertices or each other. So the 
system evolves in a 6D subspace, and we take uniform 
superpositions of identically evolving vertices/directions 
|aa), \ab), \ba), \bc), \cb), and |cc) as the basis vectors, 
e.g., 


16c) 


1 

^/M 


b£h\ue 


1 

VM- 1 


^|6^c). 

C'^b 


In terms of these basis vectors, the initial state is 


IV'o) 


1 

Vn 


M — l|aa) -|- |a6) -|- |6a) -I- a/M — l|5c) 
-h VM-l\cb) + {M- l)|cc)) , 



FIG. 2. The success probability as a function of the number 
of applications of U. The solid black line uses Uf'®’ with M = 
100, and the dashed red, dotted green, and dot-dashed blue 
curves use with M = 100, 200, and 300, respectively. 


and the quantum walk operator is 

/ cos 9 sin 0 0 0 

0 0 — cos 9 sin 9 

_ sin 0 — cos 0 0 0 

“ 0 0 0 0 

0 0 sin 0 cos 0 

V 0 0 0 0 


0 0 \ 

0 0 

0 0 

— cos 0 sin 0 ’ 

0 0 

sin 0 cos 0 / 


where cos0 = 1 — 2/M and sin0 = 2^M — 1/M. 

Note that GolV'o) = IV'o)) so to turn this quantum walk 
into a search algorithm, we use a different coin Ci for 
the marked vertices and still use Co on the unmarked 
vertices. Then the search operator is 


U = S- 


in-Y 


)Cn 




|w)(t(;| g) Cl 


Since this distinguishes the marked from the unmarked 
vertices with each application of U, it performs one oracle 
query per walk step. Two choices for Ci are common. 
The first is Cf'*’ = —Co, which causes U to become 


U = S- 





= S-{In^ Co) ■ 

■v' 

Uo 



Note Ru, = diag(—1, —1,1,1,1,1) is analogous to the 
phase flip in Grover’s algorithm nnKniEH]. The second 
choice is Cf^'^ = —Im, which gives search algorithms 
on the hypercube [3| and arbitrary dimensional periodic 
grids m- Even with these different coins for the marked 
vertices, the system still evolves in the same 6D subspace. 

Figure shows the success probability as we repeat¬ 
edly apply U with these two choices of Ci; with Cf'^, 
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the success probability stays near its initial value of 
M/N = 0(1/M), and with the success probabil¬ 

ity reaches 1/2 in 0(M) steps. To prove these behaviors, 
we simply find the eigenvectors and eigenvalues of U with 
each coin, the details of which are in Appendix With 
Cf'^, the system approximately starts in an eigenstate. 
With it approximately starts as a linear combi¬ 

nation of two eigenstates, and evolves to having equal 
probability of being in \ab) and \ba) (implying a suc¬ 
cess probability of 1/2 from the \ab) term) after ttM/2'</2 
steps. While both of these are no better than the 0(M) 
queries to guess which complete graph is marked, they 
are quadratically better than the 0 (M^) queries needed 
by the classical walk that makes one query per walk step. 


III. CONTINUOUS-TIME QUANTUM WALKS 


Now let us consider search by continuous-time quan¬ 
tum walk, which does not require the coin space, so the 
system walks on the vertices labeling an A-dimensional 
Hilbert space. The system \'ip) begins in an equal su¬ 
perposition |s^) over the vertices, and it evolves by 
Schrodinger’s equation with Hamiltonian [S] 

H =-iA-'^\w){w\, (I) 


where 7 is the jumping rate {i.e., amplitude per time), 
A is the adjacency matrix of the graph (A^ = I if j and 
j are adjacent and 0 otherwise), and w sums over the 
marked vertices. The first term effects a quantum walk 
[S] while the second term acts as an oracle [TS], with 
7 setting their relative strength. As shown in Fig. 
there are only three types of vertices. Since there are no 
directions, the system evolves in a 3D subspace spanned 
by uniform superpositions of the types of vertices: 


I a) 


1 

^/M 


H I*)- 

iGred 


\b) 


1 

^/M 


Y. i>> 

zGblue 


|c) 


I 

yjM{M - I) 


E I')- 

iGwhite 


In this 3D basis, the search Hamiltonian Q is 

(M-l+^ 1 0 \ 

H = -7 1 ^ 0 

V 0 \/M- 1 M -I ) 


Figure shows the success probability as the system 
evolves by this Hamiltonian when 7 = 1-l-l/M. As shown 
in Appendix lousing a diagrammatic approach [2^ to de¬ 
generate perturbation theory m. at this value of 7 , two 
of the eigenstates of H are approximately (|s«) ± \a))/^/2 
with eigenvalues —M— 1=F I/v/M, so the system evolves 
from |si,) to I a) in time ty/AE = tt'/M/2. 

While this appears to be a quartic speedup over the 
classical walk, it is not. The coefficient of the oracle 



FIG. 3. The success probability as a function of time for 
continuous-time search with 7 = 1-1- 1/M. The solid black, 
dashed red, and dotted green curves are M = 100, 200, and 
300, respectively. 


term in the Hamiltonian 0 is 1 , so we are counting the 
number of oracle queries and allowing an arbitrary num¬ 
ber of walk steps per oracle query. With 7 « 1, this 
corresponds to M walk steps per oracle query since the 
operator norm of A is M while the operator norm of the 
oracle term is 1. Alternatively, if we use 7 A to define 
a classical random walk, it makes each transition with 
probability per time 7 , and since there are M possible 
transitions from each vertex, the probability of making 
some transition per time is 7 M. A classical random walk 
that takes M walk steps for every oracle query optimally 
finds a marked vertex with 0 (M) queries (but 0 (M^) 
walk steps), and so the continuous-time algorithm is a 
quadratic speedup in the number of oracle queries. 


IV. MULTIPLE WALK STEPS PER ORACLE 
QUERY 

Thus we have an example where the continuous-time 
quantum walk outperforms the discrete-time quantum 
walk in search, using the oracle for Q/sfM) time rather 
than 0 (M) queries, in stark contrast to previous results 
showing that discrete-time quantum walks are faster m- 
But it does not seem fair that the discrete-time quantum 
walk only take one walk step per oracle query when the 
continuous-time algorithm is allowed to take more. To 
provide greater equity, we modify the discrete-time algo¬ 
rithm to also take multiple walk steps per oracle query 
by repeatedly applying 


U = U^R^, 

for some positive integer fc, to the initial equal superpo¬ 
sition over the coin and vertex spaces |'(/'o) = |sv) G) |sc). 
Then the number of applications of U counts the number 
of oracle queries, similar to how evolution by 0 counts 
the time that the oracle is used. A single application of 
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U acts on the initial state by 


u^RM = t/( 






M — llaa) + 


= IV'o) - -^^0 (VM -l\aa) + |o 6 )) 


y/N 

. , 2\/ M — 1 ^ 

IV'o)- ^ — Ul;\aa), 


where we used Uo\'ipo) = |V'o)- If Uq were not present, the 
second term would subtract amplitude from |aa), which 
would decrease the success probability. But we want to 
increase it instead, so we want to pick k so as to add 
success amplitude. As shown in Appendix |aa) « 
{\(j)) + \—(j)))/y/2, where \±(j)) are eigenvectors of Uq with 
eigenvalues where sint/) « y/2/M. Then 


« l^o) - + |-<^)) 

Then we pick 

^ (2n + l)7r (2n + l)7r {2n+\)'Ky/M 

(j) sin (p s/2 ’ 

with integer n, so that the exponentials equal —1, flipping 
the second term so that it adds success amplitude rather 
than decreasing it: 

« IV'o) + + 1-0)) 

, , , 2sfW^. , 

~ I0o) + —— \aa). 

With this choice of k, the algorithm takes Q{y/M) walk 
steps for each oracle query. 

As shown in Appendix O with this choice of k, two of 
the eigenvectors of 17 = U^Rw are approximately |±cr) = 
(|cc) =F \aa))/y/2 with eigenvalues where sincr « 

2/y/M. Then the initial state |0o) is approximately 

|0o) « |cc) = ^ (|cr) + |-cr)). 

Acting on this t times with UqRw, 

{U^R4 I0 o) « ^ (e^'^V) + e-^-*\-a)) . 

When 


TT TT TT y/M 

2(7 2 sin (7 4 


this becomes 

^ (i|cr) - i\-a)) 


—i\aa). 



FIG. 4. The success probability as a function of the number 
of applications of UqRw with k = [ny/M/y/2]. The solid 
black, dashed red, and dotted green curves are M = 100, 200, 
and 300, respectively. 


So the system evolves from \iPq) to the marked vertices 
with probability 1 in iry/M /4 oracle queries, as shown in 
Fig. Ill 

A classical random walk that similarly takes Q{y/M) 
walk steps for every oracle query finds a marked vertex 
with queries (but 0(M^) walk steps). Thus our 

novel algorithm achieves a cubic speedup over the corre¬ 
sponding classical random walk with the same number of 
walk steps per oracle query, and it is the first example 
of a quantum walk search with a greater-than-quadratic 
speedup over its corresponding classical random walk. 
This contrasts with more general results with quadratic 
speedups [111I23]. Of course, if we allow the classical and 
quantum walks to take different numbers of walk steps 
per oracle query, then the speedup is quadratic, so speci¬ 
fying the number of walk steps per oracle query is critical 
in comparing algorithms. 


V. CONCLUSION 


We have shown that the typical discrete-time quantum 
walk search algorithm fixes the ratio of walk steps to ora¬ 
cle queries to 1 -to-l, while the continuous-time algorithm 
allows for multiple walk steps per oracle query. When 
searching for a fully marked complete graph in the sim¬ 
plex of complete graphs, this results in the continuous¬ 
time algorithm outperforming the discrete-time’s, evolv¬ 
ing by the oracle for Q{y/M) time instead of using 0(M) 
queries. Both of these, however, are quadratic speedups 
over their corresponding classical algorithms. Providing 
greater equity, we modify the discrete-time algorithm to 
also allow multiple walk steps per oracle query. In doing 
so, the number of queries matches the Q{y/M) runtime of 
the continuous-time algorithm, and it is a cubic speedup 
over the corresponding classical random walk. 
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Appendix A: Discrete-Time Algorithm 
1. Phase Flip Coin 

With = —Co, the evolution operator U = UoRw 
is 

/— cos 0— sin0 0 0 0 0\ 

0 0 — cos 9 sin 0 0 0 

_ — sin 9 cos 0 0 0 0 0 

0 0 0 0 — cos 0 sin 0 

0 0 sin 0 cos 0 0 0 

Vo 0 0 0 sin0 COS0/ 

The eigenvalues of this are 

X, X, C. , C. , C. , C. , 

where <j)± is defined such that 

±Vm- 1 

M ’ 

- m + 1 

M ■ 

To Hnd the eigenvectors, which have the form ij) = 
(a, b, c, d, e, /)''', we solve Uip = Xtp. This yields six equa¬ 
tions: 


cos(p± = 


± sin 0 


sm 0 ± = 


\/4 — sin^ 0 


—a cos 0 — 5 sin 0 = Aa 
—CCOS0 -I- dsin0 = Xb 
—a sin 0 -I- 6 cos 9 = Xc 
—e cos 0 -I- / sin 0 = Ad 
c sin 0 -I- d cos 0 = Ae 
e sin 0 -I- / cos 0 = A/. 


Solving these yields 


, — A — COS0 — 1 —Acos0 

a = a, b= -—- a, c= -—— a, 


sin 0 


Asin 0 


, — A^ — A^ cos 0 — A cos^ 0 — cos 0 

d = -^- a, 

X sin^ 0 

—A^ cos 0 — A^ cos^ 0 — A cos 0—1 

e = --9- a, 

A2 sin^ 0 

—A^ cos 0 — A^ cos^ 0 — A cos 0—1 
A2 sin0(A — COS0) 


Let 0 = 1 . Then the (unnormalized) eigenvector corre¬ 
sponding to eigenvalue A = 1 is 

/ — 1 —COS0 — 1 —COS0 —(l-l-cos0)^ 

1 ^ B ^ B ^ 2 /I ’ 

\ sma* sm^ sin ^ 

— (1-1-cos 0 )^ —(1-1-cos 0 )^ 

sin ^0 ’ sin 0(1 — cos 0 )/ 

Plugging in for cos0 and sin0, this becomes 

V^i = ( i , -Vm^,-{m - 1 ), 

- (M- 1),-(M- 1)3/2^. 

For large M, this is dominated by the last component, 
which means that the normalized eigenvector Itpi) is ap¬ 
proximately |cc) « IV’o)- So the system approximately 
starts in an eigenstate and does not evolve significantly 
from it, as stated in the main text. 


2. SKW Coin 

With = —Im, the evolution operator U is 

/-I 0 0 0 0 0 \ 

0 0 — cos 0 sin 0 0 0 

0-10 0 0 0 

~ 0 0 0 0 — cos 0 sin 0 

0 0 sin 0 cos 0 0 0 

V 0 0 0 0 sin0 COS0/ 

The eigenvalues of this are 

_1 p-i’P- 

X, -L, C. , C. , C. , C. , 

where (j)± is defined such that 
, 1 -I- cos 0 ± a 

COS(p± = - - - 

M - 1± ^(M- l)(M-k3) 

“ TM ’ 

v^-v/5 — 2cos0 -I- cos2 0 =p a{l + cos0) 
sm(p± = - - - 

+ - 1)3/2VM -k 3 

V2M 

where 

a = \/(l -I- cos0)(5 — 3cos0) = ^a/(M — 1)(M -I- 3). 

Clearly, one of the eigenvectors corresponding to eigen¬ 
value — 1 is |aa). The remaining eigenvectors have the 
form ip = (0, a, 6 , c, d, e). To find them, we solve Utp = 
Xip, which yields five equations: 

—bcos9 -I- csin0 = Aa 
—a = Xb 

—d cos 0 -I- e sin 0 = Ac 
& sin 0 -|- c cos 9 = Xd 
d sin 0 -I- e cos 0 = Ae. 
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Solving these yields 


a = —Xb, b = b, 

, 1 — cos 9 

^ 

A sin cl 


c = 


e = 


cos 9 — X^ 
sin 6 l ’ 
1 — A^ cos 9 
A(A — COS0) 


Let 6 = 1. Then when A = the (unnormalized) 

eigenvectors are 


'4^±<j> = 


0 

-A 

1 

COS 9 —X^ 

sine 

1 —cos 9X^ 

X sin 9 
1 —cos 9X^ 

\ A(A—cos 9) f 


0 

— cos (f>^ i sin (j) 

1 

cos COS 20 -sin 20 
sin 9 “r ^ gin 9 
COS 0(1—COS 0) ■ sin 0(l+cos 0) 

sin 9 \ sin 9 

cos 20—cos 0 cos 0(1 —cos 0)—cos 0 


1+cos^ 0—2 cos 0 cos 0 




sin 20—sin 0 cos 0(l+cos 0) 
1+cos^ 0—2 cos 0 cos 0 


Then the sum and difference (times i) of the eigenvectors 
with eigenvalues and are 


/ 


' 0+0 + 0-0 = 


0 

—2 cost 

2 




*( 0+0 - 0 - 0 ) = 


2 cos 0 —cos 20 
sin 0 

2 cos 0(1 —cos 0) 
sin 0 

r) cos 20 —cos 0 cos 0(1 —cos 0)—cos 0 

V 1+cos^ 0—2 COS 0 COS 0 

0 

2 sin 0 

0 


2 sin 20 
sin 0 

2 sin 0(l+cos 0) 
sin 0 

rj sin 20—sin 0 cos 0(l+cos 0) 

A 1+cos^ 0—2 cos 0 cos 0 


With 0+, these become for large M 


■ 0 + 0 + + 0 - 0 + 


/o\ 

-2 

2 

0 

0 

Vo/ 


*( 0 + 0 + - 0 - 0 +) 


/ 0 \ 
0 
0 
0 
0 

\ 2 V 2 J 


Normalizing, we see that the initial state |0o) is approx¬ 
imately 


1 

| 0 o) « |cc) « ^(^+<^+ “ 0 - 0 +)- 

After t applications of U, the system is approximately in 
the state 


C/*l 0 o) 


i 

2 V 2 


(e*‘^+V+ 0 + - 




When 0+t = n/2, i.e., when 


TT TT TtM 

20 + 2 sin 0 + 2 ^/ 2 ’ 


then the state is 
i 


-1 


2 ^(* 0 + 0 + - (-*) 0 - 0 +) = ^( 0 + 0 + + 0 - 0 +) 


1 

7^ 


/o\ 

1 

-1 

0 

0 

V 0 


So after TrM/2^/2 applications of U, the system evolves 
from |s) to ^(|ofe) — \ba)), which has probability 1/2 of 

being at a marked vertex (from the |a 6 ) piece), as stated 
in the main text. 


Appendix B: Continuous-Time Algorithm 

The search Hamiltonian is 

/M-l + i 1 0 \ 

i? = -7 1 0 VM- 1 . 

V 0 y/W^ M -I ) 

To find the eigenvalues and eigenvectors of this, we em¬ 
ploy degenerate perturbation theory [5T]. To begin, 
we visualize the Hamiltonian diagrammatically |20j as 
a graph with three vertices, ignoring the overall factor of 
— 7 , as shown in Fig. We choose the leading-order 
Hamiltonian to exclude terms that scale less than v/M, 
so it is: 


(M+^ 0 (^\ 

= -7 0 0 y/M 

\ 0 y/M M j 


Diagrammatically, we’ve eliminated edges that scale less 


than y/M, as shown in Fig. pH From this diagram, we see 


that has three eigenvectors: one is |a), and the other 
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M-l+i 

~i 

4- 


M-1 

-<i) 

(a) 



M 



(b) 


FIG. 5. Apart from a factor of — 7 , (a) the Hamiltonian 
for continuous-time search on the simplex of complete graphs 
with one fully marked complete graph, and (b) the leading- 
order terms. 


where Har = etc. With 7 = 7 c, this 

yields 


V '/M 


Vm 

-M-1 




for large N. Solving this, the eigenstates and eigenvalues 
of H are approximately 

l« = + l«». E=-M-l-^ 

l« = ^(W-|«». E=-M-l + ^. 

Since |s„) « |r), these eigenstates are approximately 
(|s„) ± \a))/y/2, as stated in the main text. 


two are linear combinations of | 6 ) and |c). Specifically, 

they are Appendix C: Multiple Walk Steps per Oracle Query 


i), E = —'y i M — 


-VM-^ M - ^M(M + 4) 

- 2 -1^) + E = -7_™-X™™__ 

-TM -b VM + 4.^. ^ M -h JM{M + 4) 

---p) + |c), E = --f -- 


Note the third state is approximately 


1 

Vm 


+ |c) 


h T. T. !■> 

iGblue ^MGwhite 

E I*) = 1 ^)- 

i^red 


where |r) is the uniform superposition over the unmarked 
vertices. Setting the first and third eigenvalues equal (so 
that |a) and |r) are degenerate), we get the critical 7 : 

~ -M + ^M{M + A)^ "^M' 

The perturbation restores terms of 0(1), which in 
the diagram includes the edge of weight 1 between |a) and 
\b), the —1 in |a)’s self-loop, and the —1 in |c)’s self-loop. 
Away from the critical 7 , the perturbation does not 
change the eigenvectors significantly, and so the initial 
state is approximately an eigenstate and fails to evolve 
beyond a global, unobservable phase. Near 7 c, however, 
the perturbation causes two linear combinations of |a) 
and |r). 


\tp) = aa\a) + ar\r), 

to be eigenstates of the perturbed system. The coeffi¬ 
cients can be found by solving 

f^aa ^ar\ f^a\ 

HrrJ \arj ~ ^ \ar J ’ 


Let us start by finding the eigenvalues and eigenstates 
of the quantum walk operator 


/cosd 

sind 

0 

0 

0 


0 

0 

— cosd 

sind 

0 

0 

sin d 

— cosd 

0 

0 

0 

0 

0 

0 

0 

0 

— cosd 

sind 

0 

0 

sin d 

cosd 

0 

0 

V 0 

0 

0 

0 

sind 

cos 0/ 


This has eigenvalues 

l,-l,e*‘^+,e-*'^+,e*'^-,e-*‘^-, 

where 

±l-l-cos0 -^(1 =F cos 9){3 ± cos9) 

cos 0 ± =- - -, sm 0 ± =- - -. 

To find the eigenstates, which have the form ip = 
(a, b, c, d, e, /)''', we solve Uip = Xip. This yields six equa¬ 
tions: 


a cos 0 + b sin 0 = Xa 
—CCOS0 + dsm0 = Xb 
a sin 0 — b cos 0 = Xc 
—e cos 0 + f sm0 = Xd 
c sin d -I- d cos 0 = Xe 
e sin 6 * -I- / cos 0 = Xf. 


Solving these yields 


, A — cos 0 1 — A cos 0 

a = a, b = -^—-—a, c = —r—^— - —a. 


sin 0 


Asind 


, A^ — A^ cos 0 — X cos^ 0 + cos 0 

d = -^- a, 

X sin^ 0 

A^ cos 0 — X^ cos^ 0 — X cos d -I-1 

e =-^T 

A^ sin^ 0 

A^ cos d — A^ cos^ d — A cos d -I-1 
A^ sind(A — cosd) 


































Let a = 1. Plugging in for A, the (unnormalized) eigen¬ 
vectors with eigenvalues 1 and —1 are; 


V'l 

V’-I = 




1 , -Vm - 1 , -Vm -1, 1 , 1 , 




-1 


VM-1 


1 


When A = the (unnormalized) eigenvectors are 


i’±4> 


cos 0—cos t 


2 

cos 20—cos 0 


± i 


zb i 

Ti . 

sin(20) —2 sin 0 cos 9 


sin 9 
sin 0 


sin^9 sin^9 

cos20—cos^0 _ • sin(20) —2 sin 0 cos 0 
sin^ 9 ^ 3 ^ 

(1 — 2 cos 20) cos 0—cos 0 cos^ 0+cos 30+cos'^ 9 ■ sin(30) —2 sin(20) cos 0+sin 0 cos^ 9 

sin 0(1 —2 cos 0 cos 0+cos2 0) sin 0(1 —2 cos 0 cos 0+cos2 0) > 


Plugging in for (/>+ and </)_, as well as for cos^ and sin^, 
this yields eigenvectors 




'4’±<p+ 


( 2 ^_\ 

1 _|- W 2 M -1 

VM-1 VM-l 
1 1 

VM-l ' VM-1 

2 ^-VlMEI 

M-1 ^ ‘ M-1 
-1 ■ V2M-1 

M-1 “1“ ‘ M-1 

\ VM-1 / 


l±<^-) 


Then 


and 


/I 1 ± 1 * 

Vv^’ 20 W 2 ’ 20 W^ 2 ’ 

-1 ^ i -1 i -1 y 

2y/2M 2v^’ 2v^M ^ 2v^’ ^/W/ 

/ 1 -1 ^ i -1 i 

\M'2^/M 2v^’2v{M ^ 2v^’ 

—1 i —I i ^ y 


|aa) « ^ ( 10 +) + 1 - 0 +)), 


'4’±<P- 


( 


1 

2 


V 


2 \ 

-VM - 1 ± WM + 1 

(M - 1) =F VM2 - 1 
(M - 1) ± VM2 - 1 
2 VM- 1 / 


as stated in the main text. 

Now let us find the eigenvalues and eigenvectors of 
U^Rw when k = {2n + /\/2 (which we assume 

from now on). Using the above eigenstates, we can de¬ 
fine the matrices 

P=(|l) 1-1) 10+) 1-0+) |0_) |-0_)) 


Let us take large M so that — 1« M, 2M — 

1 « 2M, and — 1 « M^. Then all the normalized 
eigenstates are approximately 


and 

D'^ = diag ( 1 , (-l)^ - 1 , - 1 , 6 *'='^-, 6 -*'='^-) 



1111 1 

0 m’ M’ M’ 0 m’ 0 m 



T 


1 - 1 ) 


1-1-11 1 -1 y 

0W’ 02’ 02’ 0W’ 0W’ 02M/ 


that diagonalize C/q. Then 

is, up to terms O( 1 / 0 M), 
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1 

l+(-l)'' 

-l-(-l)'' 

0 

0 

2 

2VM 

2VM 


! + (-!)'= 

Hi -(-in 


— l—(—1)^+2 sin(fc(/)_) 

l-(-l)*“+2cos(fc0_) 

0 

2VM 

2VM 

2VM 

! + (-!)'= 

n-i-(-in 

1 (-l + (-l)fc) 

l-{-l)’‘+2cos{k(p-) 

— 1—(—1)^—2 sin(fe^_) 

0 

2VM 

2VM 

2VM 

0 

l+( — 1)^+2 sin(fc0_) 

l-(-l)'“+2cos(fe0_) 

cos{k4>-) 

— sin(fc(/)_) 

1 —cos(fc0_ )+sin(/i:(/i_ ) 

2^/M 

2VM 


0 

-l-|-(-l)''-2cos(fc0_) 

— 1 —( —1)^+2 sin(k(p—) 

sin{k(j)-) 

cos{k(j)-) 

l—cos{k4>— ) — sm{k(l )—) 

2VM 

2^/M 

Vm 


0 

0 

l — cos{k(p—)—sin{k(f )-) 

l — cos{k(f)- )-\-sin{k(p—) 

1 


vVtf 



To find the eigenvalues and eigenvectors of this, we use 
degenerate perturbation theory. The leading-order terms 
are 


{U^Rn,Y°^ 


/I 
0 

0 i (-!-(- 


0 

0 

Vo 


5(1-(-!)") 5(1+ (-!)") 


0 0 0 \ 

0 0 0 

0 0 0 

cos{k(j)-) — sin(fc(/)_) 0 
sin(A:^_) cos{k(j)-) 0 

0 01/ 


r 


Let us separately consider the cases when k is even and 
when k is odd. 

When k is even, the leading-order terms of U^Rw are 


[U^R 


0 )k even 




/I 0 0 0 0 

0010 0 

0-100 0 0 

000 cos{k(j>-) —sm{k(j>-) 0 
0 0 0 sin(fc(/)_) cos{k4>-) 0 

\0 0 0 0 0 ij 


which has eigenvectors and eigenvalues 
^ {-i\ab) + \ba )), i 


V2 


{i\ab) + |6a)), 


|aa), 

|cc). 


^ (-i|6c) -h \cb )), 
^ (z|6c) -h \cb )), 


^—ik4>— 




The eigenstates we care about are |aa) and |cc), which are 
degenerate. The perturbation {UqRw)^^'> restores terms 
of 0(1/-\/M) in the matrix, and two linear combinations 
of |aa) and |cc), 

aa\aa) + adcc). 


are eigenstates of the perturbed system. The coefficients 
can be found by solving 


Hrr. Hr. 


= E 


a. 


where = (aa|(Cfo+ {UqRw)^^'^\cc), etc. This 
yields 



= E 


Solving this, the eigenvectors and eigenvalues of C/q 
for even k are 


|±cr) = ^ (T|aa) + |cc)), E± = l± -^i « 

where sina « 2/'/M, as stated in the main text. 
What about when k is odd? We get 




(0) 

0 odd 


{U^R 


/I 0 0 0 0 0\ 

0 1 0 0 0 0 

00-10 0 0 
0 0 0 cos{k4>-) —sin{k4>-) 0 

0 0 0 sm{k(j)-) cos{k(j)-) 0 

Vo 0 0 0 0 1/ 
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which has eigenvectors and eigenvalues 

\ba), -1 

|aa), 1 

\ab), 1 

|cc), 1 

^{-^\bc) + \cb)), 

^ (*|6c) + |c6)), 


where Hab = (aa|(C/Q + (C/getc. This 

yields 

0 1 0 = S at 

7F ^ / \'^=/ 

Solving this, the eigenvectors are 
\ab), 1 

and 



The three eigenstates we care about are the degenerate 
states |aa), \ab), and |cc). With the 0(1/-\/M) pertur¬ 
bation {UQRyj)^^\ three superpositions of these states 
ttalaa) -I- ab\ab) + adcc) become eigenstates of the per¬ 
turbed system, where the coefficients can be found by 
solving 

( R-aa R-ah Rac\ ( rra\ ( rra\ 

Rba Rbb Rbc 1 I 1 — I ’ 

Rea Rcb Rcc j j J 


|±cr) = ^ (T|aa)-f |cc)), 1 ± 

where sin cr « 2I'/M. So we get the same result as for 
even k. While not relevant for our results, note that the 
term \ab) can have an affect after the first peak in success 
probability, causing a double peak, which the analytics 
reveal when higher-order terms are included in the cal¬ 
culation. 
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